A 1.1 Mb duplication CNV on chromosome 17 contributes to skeletal muscle development in Boer goats.

The Boer goat is one of the top meat breeds in modern animal husbandry and has attracted widespread attention for its unique growth performance. However, the genetic basis of muscle development in the Boer goat remains obscure. In this study, we identified specific structural variants in the Boer goat based on genome-wide selection signals and analyzed the basis of the molecular heredity of related candidate genes in muscle development. A total of 9 959 autosomal copy number variations (CNVs) were identified through selection signal analysis in 127 goat genomes. Specifically, we confirmed that the highest signal CNV (HSV) was a chromosomal arrangement containing an approximately 1.11 Mb (CHIR17: 60062304-61171840 bp) duplicated fragment inserted in reverse orientation and a 5 362 bp deleted region (CHIR17:60145940-60151302 bp) with overlapping genes (e.g., ARHGAP10, NR3C2, EDNRA, PRMT9, and TMEM184C). The homozygous duplicated HSV genotype (+/+) was found in 96% of Boer goats but was not detected in Eurasian goats and was only detected in 4% of indigenous African goats. The expression network of three candidate genes ( ARHGAP10, NR3C2, and EDNRA) regulating dose transcription was constructed by RNA sequencing. Results indicated that these genes were involved in the proliferation and differentiation of skeletal muscle satellite cells (SMSCs) and their overexpression significantly increased the expression of SAA3. The HSV of the Boer goat contributed to superior skeletal muscle growth via the dose effects of overlapping genes.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene.

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The FST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

Skin transcriptome analysis identifies the key genes underlying fur development in Chinese Tan sheep in the birth and Er-mao periods.

Hair follicle development in Tan sheep differs significantly between the birth and Er-mao periods, but the underlying molecular mechanism is still unclear. We profiled the skin transcriptomes of Tan sheep in the birth and Er-mao periods via RNA-seq technology. The Tan sheep examined consisted of three sheep in the birth period and three sheep in the Er-mao period. A total of 364 differentially expressed genes (DEGs) in the skin of Tan sheep between the birth period and the Er-mao period were identified, among which 168 were upregulated and 196 were downregulated. Interestingly, the FOS proto-oncogene (FOS) (fold change = 22.67, P value = 2.15*10^-44) was the most significantly differentially expressed gene. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the FOS gene was significantly enriched in the signaling pathway related to hair follicle development. Immunohistochemical analysis showed that the FOS gene was expressed in the skin of Chinese Tan sheep at the birth and Er-mao periods, with significantly higher expression in the Er-mao period. Our findings suggest that the FOS gene promotes hair follicle development in Tan sheep.

Progress on genetic mapping and genetic mechanism of cattle and sheep horns.

Horns are cranial appendages, which are unique in ruminants. Cattle (Bos taurus) and sheep (Ovis aries) cranial appendages exhibit various forms of morphology, including wild-type two-horn phenotype, polled phenotype and scur phenotype. These animals provide an ideal model for studies on the underlying relationship between quality and quantitative traits of cattle and sheep horn and the molecular mechanisms of horn phenotype as a polygenic regulation for quality traits. In recent years, some research progresses of cattle and sheep horns are successively reported, which helps us better understand the evolutionary origin of new organ, the effects of natural selection, sex selection and artificial selection on horn phenotypes. In this review, we introduce in details the recent advances on the research of horn traits in cattle and sheep, and summarize the genetic mapping of multi-horned phenotypes, the genetic mapping of polled locus, and studies on scur phenotype. Moreover, we discuss potential problems in such research, thereby providing a reference for investigation on the genetic mechanisms of horn traits in ruminants.

Evolutionary relationship and population structure of domestic Bovidae animals based on MHC-linked and neutral autosomal microsatellite markers.

Major histocompatibility complex (MHC) genes are critical for disease resistance or susceptibility responsible for host-pathogen interactions determined mainly by extensive polymorphisms in the MHC genes. Here, we examined the diversity and phylogenetic pattern of MHC haplotypes reconstructed using three MHC-linked microsatellite markers in 55 populations of five Bovidae species and compared them with those based on neutral autosomal microsatellite markers (NAMs). Three-hundred-and-forty MHC haplotypes were identified in 1453 Bovidae individuals, suggesting significantly higher polymorphism and heterozygosity compared with those based on NAMs. The ambitious boundaries in population differentiation (phylogenetic network, pairwise FST and STRUCTURE analyses) within and between species assessed using the MHC haplotypes were different from those revealed by NAMs associated closely with speciation, geographical distribution, domestication and management histories. In addition, the mean FST was significantly correlated negatively with the number of observed alleles (NA), observed (HO) and expected (HE) heterozygosity and polymorphism information content (PIC) (P < 0.05) in the MHC haplotype dataset while there was no correction of the mean FST estimates (P> 0.05) between the MHC haplotype and NAMs datasets. Analysis of molecular variance (AMOVA) revealed a lower percentage of total variance (PTV) between species/groups based on the MHC-linked microsatellites than NAMs. Therefore, it was inferred that individuals within populations accumulated as many MHC variants as possible to increase their heterozygosity and thus the survival rate of their affiliated populations and species, which eventually reduced population differentiation and thereby complicated their classification and phylogenetic relationship inference. In summary, host-pathogen coevolution and heterozygote advantage, rather than demographic history, act as key driving forces shaping the MHC diversity within the populations and determining the interspecific MHC diversity.

Convergent genomic signatures of high-altitude adaptation among domestic mammals.

Abundant and diverse domestic mammals living on the Tibetan Plateau provide useful materials for investigating adaptive evolution and genetic convergence. Here, we used 327 genomes from horses, sheep, goats, cattle, pigs and dogs living at both high and low altitudes, including 73 genomes generated for this study, to disentangle the genetic mechanisms underlying local adaptation of domestic mammals. Although molecular convergence is comparatively rare at the DNA sequence level, we found convergent signature of positive selection at the gene level, particularly the EPAS1 gene in these Tibetan domestic mammals. We also reported a potential function in response to hypoxia for the gene C10orf67, which underwent positive selection in three of the domestic mammals. Our data provide an insight into adaptive evolution of high-altitude domestic mammals, and should facilitate the search for additional novel genes involved in the hypoxia response pathway.

Identification of a Goat Intersexuality-Associated Novel Variant Through Genome-Wide Resequencing and Hi-C.

Background: Polled intersex syndrome (PIS) leads to reproductive disorders in goats and exerts a heavy influence on goat breeding. Since 2001, the core variant of an 11.7 kb deletion at ~129 Mb on chromosome 1 (CHI1) has been widely used as a genetic diagnostic criterion. In 2020, a ~0.48 Mb insertion within the PIS deletion was identified by sequencing in XX intersex goats. However, the suitability of this variation for the diagnosis of intersex goats worldwide and its further molecular genetic mechanism need to be clarified. Results: The whole-genome selective sweep of intersex goats from China was performed with whole-genome next-generation sequencing technology for large sample populations and a case-control study on interbreeds. A series of candidate genes related to the goat intersexuality phenotype were found. We further confirmed that a ~0.48 Mb duplicated fragment (including ERG and KCNJ15) downstream of the ~20 Mb PIS region was reversely inserted into the PIS locus in intersex Chinese goats and was consistent with that in European Saanen and Valais black-necked goats. High-throughput chromosome conformation capture (Hi-C) technology was then used to compare the 3D structures of the PIS variant neighborhood in CHI1 between intersex and non-intersex goats. A newly found structure was validated as an intrachromosomal rearrangement. This inserted duplication changed the original spatial structure of goat CHI1 and caused the appearance of several specific loop structures in the adjacent ~20 kb downstream region of FOXL2. Conclusions: Results suggested that the novel complex PIS variant genome was sufficient as a broad-spectrum clinical diagnostic marker of XX intersexuality in goats from Europe and China. A series of private dense loop structures caused by segment insertion into the PIS deletion might affect the expression of FOXL2 or other neighboring novel candidate genes. However, these structures require further in-depth molecular biological experimental verification. In general, this study provided new insights for future research on the molecular genetic mechanism underlying female-to-male sex reversal in goats.

Identification of lncRNAs by RNA Sequencing Analysis During in Vivo Pre-Implantation Developmental Transformation in the Goat.

Pre-implantation development is a dynamic, complex and precisely regulated process that is critical for mammalian development. There is currently no description of the role of the long noncoding RNAs (lncRNAs) during the pre-implantation stages in the goat. The in vivo transcriptomes of oocytes (n = 3) and pre-implantation stages (n=19) at seven developmental stages in the goat were analyzed by RNA sequencing (RNA-Seq). The major zygotic gene activation (ZGA) event was found to occur between the 8- and 16-cell stages in the pre-implantation stages. We identified 5,160 differentially expressed lncRNAs (DELs) in developmental stage comparisons and functional analyses of the major and minor ZGAs. Fourteen lncRNA modules were found corresponding to specific pre-implantation developmental stages by weighted gene co-expression network analysis (WGCNA). A comprehensive analysis of the lncRNAs at each developmental transition of high correlation modules was done. We also identified lncRNA-mRNA networks and hub-lncRNAs for the high correlation modules at each stage. The extensive association of lncRNA target genes with other embryonic genes suggests an important regulatory role for lncRNAs in embryonic development. These data will facilitate further exploration of the role of lncRNAs in the developmental transformation in the pre-implantation stage.

Genome-wide selection signatures analysis of litter size in Dazu black goats using single-nucleotide polymorphism.

Litter size is considered to be the most important index for estimating domestic animal productivity. The number of indigenous goats in China with higher litter sizes than those of commercial breeds in other countries may be helpful for accelerating genetic improvements in goat breeding. We performed a genome-wide selective sweep analysis of 31 Dazu black goats with extreme standard deviation in litter size within the third fetus to identify significant genomic regions and candidate genes through different analyses. The analysis identified a total of 33,917,703 variants, including 32,262,179 SNPs and 1,655,524 indels. In addition, two novel candidate genes (LRP1B and GLRB), which are related to litter size, were obtained with π, Tajima's D, πA/πB, and F ST at the individual level with a 95% threshold for each parameter. These two genes were annotated in five GO terms (localization, binding, macromolecular complex, membrane part, and membrane) and two pathways (long-term depression and neuroactive ligand-receptor interaction pathway). Regarding the result of linkage disequilibrium (LD) analysis, in LRP1B and GRID2, the high-yield Dazu black goats exhibit significantly different LD patterns from low-yield goats. Litter size variability has low heritability and is related to multiple complex factors found in domestic animals. Obtaining a clear explanation and significant signal by genome-wide selective sweep analysis with a small sample size is difficult. However, we investigated some candidate genes, particularly LRP1B and GLRB, which may provide useful information for further research.

Genetic diversity estimation of Yunnan indigenous goat breeds using microsatellite markers.

BACKGROUND: To assess the genetic diversity of seven Yunnan indigenous goat populations (Fengqing hornless goat, Mile red-bone goat, Longling goat, Ninglang black goat, Black-bone goat, Yunling black goat, and Zhaotong goat), their population structures were investigated using 20 microsatellite markers. RESULTS: The results indicated that the genetic diversity of these goats was rich. The observed heterozygosity ranged from 0.4667 ± 0.0243 to 0.5793 ± 0.0230, and the mean number of alleles ranged from 4.80 ± 1.61 and 4.80 ± 1.64 to 6.20 ± 2.93. The population structure analysis showed that these seven goat populations were separated into two clusters, consistent with the results from phylogenetic networks, pairwise differences, and STRUCTURE analyses. We speculate that this may have been caused by natural geographical isolation, human migration and economic and cultural exchanges. We suggest removing CSRD247 and ILSTS005, two loci identified to be under positive selection in the present study, from the microsatellite evaluation system of goats. CONCLUSIONS: The present study may provide a scientific basis for the conservation and utilization of Yunnan indigenous goats.

Current genetic diversity in eight local Chinese sheep populations.

China has numerous local domestic sheep breeds. In this study, the genetic diversity of eight sheep populations was estimated using 17 microsatellites. Knowledge of such diversity provides novel insight into the degree of breed protection needed and the prediction of hybrid advantage. In total, 17 microsatellites were genotyped in 186 individuals from eight populations. The mean number of alleles (± SD) ranged from 3.71 ± 1.36 in Zhaotong sheep to 11.94 ± 3.58 in small-tailed Han sheep. The observed heterozygote frequency (± SD) within a population ranged from 0.482 ± 0.025 in Zhaotong sheep to 0.664 ± 0.023 in Tibetan sheep. In addition, using pairwise difference (FST) analysis, the highest within-population diversity was observed in Tibetan sheep (πX = 12.8098) and small-tailed Han (πX = 12.67873), and the lowest diversity was observed in Zhaotong sheep (πX = 7.90337). The results for genetic divergence between populations indicated that the populations were significantly different (P < 0.05) based on the average number of pairwise differences between populations (πXY) and the corrected average pairwise differences. Both phylogenetic networks and structure analysis showed that these eight populations were separated into three clusters in accordance with their geographical habitat, except Tibetan and Hu sheep. In short, we genotyped eight local Chinese sheep populations using 17 microsatellites, and the results indicated that their current genetic diversity is decreasing and that new conservation strategies are needed. In addition, significant genetic differences between populations could be used in cross breeding.

Genetic diversity of the Chinese goat in the littoral zone of the Yangtze River as assessed by microsatellite and mtDNA.

The objective of this study was to assess the genetic diversity and population structure of goats in the Yangtze River region using microsatellite and mtDNA to better understand the current status of those goat genetic diversity and the effects of natural landscape in fashion of domestic animal genetic diversity. The genetic variability of 16 goat populations in the littoral zone of the Yangtze River was estimated using 21 autosomal microsatellites, which revealed high diversity and genetic population clustering with a dispersed geographical distribution. A phylogenetic analysis of the mitochondrial D-loop region (482 bp) was conducted in 494 goats from the Yangtze River region. In total, 117 SNPs were reconstructed, and 173 haplotypes were identified, 94.5% of which belonged to lineages A and B. Lineages C, D, and G had lower frequencies (5.2%), and lineage F haplotypes were undetected. Several high-frequency haplotypes were shared by different ecogeographically distributed populations, and the close phylogenetic relationships among certain low-frequency haplotypes indicated the historical exchange of genetic material among these populations. In particular, the lineage G haplotype suggests that some west Asian goat genetic material may have been transferred to China via Muslim migration.

'Ectopic' suprasellar type IIa PRL-secreting pituitary adenoma.

BACKGROUND: Ectopic pituitary adenomas (EPAs) are rare, and the suprasellar cistern seems to be the most common location. At this time, no detailed original classification, diagnosis, or treatment protocols for suprasellar pituitary adenomas (SPAs) have been described. CASE DESCRIPTION: A 19-year-old man showed visual disturbances and lack of libido for 3 years, he suffered a sharp decline in vision with only light perception in the last week. Magnetic resonance imaging scans revealed a large suprasellar cystic lesion with a normal pituitary in the sella turcica. Endocrinological findings showed an extremely high prolactin level of 1250 ng/mL. Because of the sharp decline in vision, the patient underwent total removal of the suprasellar lesion using a transfrontal interhemispheric approach. The tumor pedicle originated in the lower pituitary stalk without any connection to the anterior pituitary gland in the sella turcica, while the diaphragma sellae was incomplete. Clinical and endocrinological cure criteria were fulfilled and postoperative pathology confirmed a prolactin-secreting pituitary adenoma. CONCLUSION: Ectopic suprasellar pituitary adenomas (ESPAs) are extremely rare intracranial extracerebral tumors. SPAs can be classified into three types according to their origin and their relationship with surrounding tissue. Only type III is theoretically a true ectopic, based on previous reports. Thus, ESPAs are uncommon compared to other EPAs. Our case is the first reported case of a type IIa 'E'SPA and the first description of this subtype classification until now. The pars tuberalis may be different from the pars distalis, and each subtype of adenohypophyseal cells may have different migration characteristics, which leads to different proportions of each hormone-secreting subtype in SPAs and EPAs. Transsphenoidal surgery is minimally invasive, but transcranial surgery may remain a universal option for the treatment of suprasellar lesions.

Downregulation of miR-199b promotes the acute spinal cord injury through IKKß-NF-κB signaling pathway activating microglial cells.

Inflammatory response played an important role in the progression of spinal cord injury (SCI). Several miRNAs were associated with the pathology of SCI. However, the molecular mechanism of miRNA involving in inflammatory response in acute SCI (ASCI) was poorly understood. Sprague-Dawley (SD) rats were divided into 2 groups: control group (n=6) and acute SCI (ASCI) group (n=6). The expression of miR-199b and IκB kinase ß-nuclear factor-kappa B (IKKß-NF-κB) signaling pathway were evaluated by quantitative reverse transcription-PCR (qRT-PCR) in rats with ASCI and in primary microglia activated by lipopolysaccharide (LPS). We found that downregulation of miR-199b and activation of IKKß/NF-κB were observed in rats after ASCI and in activated microglia. miR-199b negatively regulated IKKß by targeting its 3'- untranslated regions (UTR) through using luciferase reporter assay. Overexpression of miR-199b reversed the up-regulation of IKKß, p-p65, tumor necrosis factor-α (TNF-α) and interleukin-1ß (IL-1ß) in LPS-treated BV2 cells assessed by western blotting analysis. In addition, BMS-345541 reversed the up-regulation effects of miR-199b inhibitor on the expression of TNF-α and IL-1ß. In the SCI rats, overexpression of miR-199b attenuated ASCI and decreased the expression of IKKß-NF-κB signaling pathway and TNF-α and IL-1ß. These results indicated that miR-199b attenuated ASCI at least partly through IKKß-NF-κB signaling pathway and affecting the function of microglia. Our findings suggest that miR-199b may be employed as therapeutic for spinal cord injury.

Expression profiling and functional characterization of miR-192 throughout sheep skeletal muscle development.

MicroRNAs (miRNAs) are evolutionarily conserved, small, non-coding RNAs that have emerged as key regulators of myogenesis. Here, we examined the miRNA expression profiles of developing sheep skeletal muscle using a deep sequencing approach. We detected 2,396 miRNAs in the sheep skeletal muscle tissues. Of these, miR-192 was found to be up-regulated in prenatal skeletal muscle, but was down-regulated postnatally. MiR-192 expression also decreased during the myogenic differentiation of sheep satellite cells (SCs). MiR-192 overexpression significantly attenuated SCs myogenic differentiation but promoted SCs proliferation, whereas miR-192 inhibition enhanced SCs differentiation but suppressed SCs proliferation. We found that miR-192 targeted retinoblastoma 1 (RB1), a known regulator of myogenesis. Furthermore, knockdown of RB1 in cultured cells significantly inhibited SCs myogenic differentiation but accelerated SCs proliferation, confirming the role of RB1 in myogenesis. Taken together, our findings enrich the ovine miRNA database, and outline the miRNA transcriptome of sheep during skeletal muscle development. Moreover, we show that miR-192 affects SCs proliferation and myogenic differentiation via down-regulation of RB1.

Conservation genetics in Chinese sheep: diversity of fourteen indigenous sheep (Ovis aries) using microsatellite markers.

The domestic sheep (Ovis aries) has been an economically and culturally important farm animal species since its domestication around the world. A wide array of sheep breeds with abundant phenotypic diversity exists including domestication and selection as well as the indigenous breeds may harbor specific features as a result of adaptation to their environment. The objective of this study was to investigate the population structure of indigenous sheep in a large geographic location of the Chinese mainland. Six microsatellites were genotyped for 611 individuals from 14 populations. The mean number of alleles (±SD) ranged from 7.00 ± 3.69 in Gangba sheep to 10.50 ± 4.23 in Tibetan sheep. The observed heterozygote frequency (±SD) within a population ranged from 0.58 ± 0.03 in Gangba sheep to 0.71 ± 0.03 in Zazakh sheep and Minxian black fur sheep. In addition, there was a low pairwise difference among the Minxian black fur sheep, Mongolian sheep, Gansu alpine merino, and Lanzhou fat-tailed sheep. Bayesian analysis with the program STRUCTURE showed support for 3 clusters, revealing a vague genetic clustering pattern with geographic location. The results of the current study inferred high genetic diversity within these native sheep in the Chinese mainland.

High occurrence of mitochondrial heteroplasmy in nepalese indigenous sheep (Ovis aries) compared to chinese sheep.

Heteroplasmy due to length polymorphism with tandem repeats in mtDNAs within individual was hardly studied in domestic animals. In the present study, we identified intra-individual length variation in the control region of mtDNAs in Nepalese sheep by molecular cloning and sequencing techniques. We observed one to four tandem repeats of a 75-bp nucleotide sequences in the mtDNA control region in 45% of the total Nepalese sheep sampled in contrast to the Chinese sheep, indicating that the heteroplasmy is specific to Nepalese sheep. The high rate of heteroplasmy in Nepalese sheep could be a resultant of the mtDNA mutation and independent segregation at intra-individual level or a strand slippage and mispairing during the replication.

A complete mitochondrial genome of Dazu Black goat.

Dazu Black goat is an indigenous goat genetic resource in Southwest of China. Here, we describe its complete mitochondrial genome sequence. The mitogenome is 16,641 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.5%, T: 27.3%, C: 26.1% and G: 13.1%. The complete mitogenome of the indigenous goat could provide important data to further explore the taxonomic status of the subspecies and also provide a starting point for further phylogenetic studies.

Complete mitochondrial genome of Chuanzhong black goat in southwest of China (Capra hircus).

The Chuanzhong black goat (Capra hircus) is a breed native to southwest of China. Its complete mitochondrial genome is 16,641 nt in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.5%, T: 27.3%, C: 26.1%, and G: 13.1%. The complete mitogenome of the Chinese indigenous breed of goat could provide a basic data for further phylogenetics analysis.

Characteristics of the mitochondrial genome of four native goats in China (Capra hircus).

Here, we describe the complete mitochondrial genome sequences of Jining Gray goat, Fushun black goat, Youzhou black-skin goat, and Hechuan white goat. The mitogenome of those four goats consisted of 16,640 nt, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. The complete mitogenome of these four local breeds of Chinese native goats could provide an important data to further breed improvement and animal genetics resource conservation in China.